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ABSTRACT 

For the 2dFGRS we study the properties of voids and of fainter galaxies within voids 
that are defined by brighter galaxies. Our results are compared with simulated galaxy 
catalogues from the Millenium simulation coupled with a semianalytical galaxy for- 
mation recipe. We derive the void size distribution and discuss its dependence on 
the faint magnitude limit of the galaxies defining the voids. While voids among faint 
galaxies are typically smaller than those among bright galaxies, the ratio of the void 
sizes to the mean galaxy separation reaches larger values. This is well reproduced in 
the mock galaxy samples studied. We provide analytic fitting functions for the void 
size distribution. Furthermore, we study the galaxy population inside voids defined 
by objects with Bj — 51og/i < —20 and diameter larger than 10 h"-^ Mpc . We find 
a clear bimodality of the void galaxies similar to the average comparison sample. We 
confirm the enhanced abundance of galaxies in the blue cloud and a depression of the 
number of red sequence galaxies. There is an indication of a slight blue shift of the 
blue cloud. Furthermore, we find that galaxies in void centers have higher specific star 
formation rates as measured by the r] parameter. We determine the radial distribution 
of the ratio of early and late type galaxies through the voids. We find and discuss 
some diff'erences between observations and the Millenium catalogues. 
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The large-scale galaxy distribution is highly inhomoge- 
neous. We observe groups, clusters and superclusters of 
galaxies and large voids. During last decades, much at- 
tention was paid on the analysis of bound structures as 
groups and clusters. Recently, new superclusters catalogues 
were constructed from the 2d FGRS and compared with 
large cosmological simulations (|Einasto et al. I l2007bm . In 
a complement, there are large regions in the universe 
without bright galaxies, so called cosmic voids. Early on 
ve ry large voids over 50 Mpc diameter were foun d 
bv lGregorv fc Thompson! l| 19781 ) and lKirshner et al.l l|l98lh . 
More common are voids with diameters of about 10 Mpc 
that fill most of cosmic space. The explanation of such struc- 
tures is not obvious. According to the standard paradigm of 
cosmological structure formation, negative potential wells 
from primordial inhomogeneities attract all matter in bound 
structures. In the same way, positive potential perturbations 
expel matter, but observed voids are too large for completely 
emptying. Therefore, in addition to the dilution of mat- 
ter, the galaxy formation probab ility should be suppresse d 
in underdense regio ns, cp. e.g. Lee fc ShandarinI ll 19981): 
iMadsen et al] ll 19981). Recently , [Sheth fc van de WeygaertI 
1 2004 ) and iFurlanetto fc PiranI (|2006l ) applied these ideas 



within the excursion set formalism of gravitational insta- 
bility. These analytical theories derived void size distri- 
butions that are peaked typically at diameters below 10 
Mpc which seem t o be smaller than observed void 
sizes, cp. iMiiller et al.1 (|2000| ). and void sizes in CDM- 
simulati ons coupled with se mianalytical galaxy formation 
models. iBenson et al.1 l|2003l ). 

Voids were routinely ide ntified in all wide-field 
redshift surveys as the CfA (Ide Lapparent et al.l Il986l : 



IVogeley et al 



1994), the SSRS2 ([El-Ad 



llHoyle fc Vogeleyl 12004 



m - 

the LCRS (iMiiller et al] l2000l: lArbabi-Bidgoli fc MiiUe 



PiraiT ''199 



|2002|). the IRAS-survey fEl-Ad fc Pira^i200a) the 2dFGRS 



Croton et al 



20041: iPatiri et al] 



l2006l ). the SPSS (|Roias et al.ll2004 120051 : iPatiri et al.ll2006l ). 
and the DEEP 2 survey with an a nalysis of voids up to 
redshift z ~ 1 (|Conrov et al.ll2005l ). However, many void 
searches are only devoted to the identification of large voids, 
other void finders depend on special procedures as firstly 
identifying wall galaxies by an overdensity cri terion and then 
looking fo r voids bo unded by wall galaxies (lEl-Ad fc PiranI 
,1997; Ho vle fc Vogeley 2004). Furthermore, the void search 
depends on the galaxy sample used for defining voids, in 
particular on the limiti ng magni t ude o f the galaxy sample. 
In an influential paper, iPeeblesI l|200ll ) derived from near- 
est neighbor statistics that galaxies of different brightness 
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respect the same voids. He claimed that this contradicts 
the standard CDM scenario of galaxy and structure for- 
mation that seem to predict a hierarchy of galactic struc- 
tures with smaller structure for fainter objects sitting in less 
massive dark matter halos, i.e also smaller voids for fainter 
object s. In a follow up theoretical study, iMathis fc W"hit3 
l|2002l ) showed from high-resolution simulation that voids 
defined by bright galaxies are also underdense in faint galax- 
ies, i.e. that bright and faint galaxies respect similar voids. 
We want to take up this question since our e arlier studies 
of vo i ds in LCRS and in LCD M-simulations (|Miiller et al.l 
I2OOOI : lArbabi-Bidgoh fc Mulleij [2b02') showed a dependence 
of the void size distribution on the brightness limit of the 
galaxies under study and a characteristic void size scaling 
relation. iBenson et al.l (|2003t ) confirmed this scaling relation 
in simulated galaxy distributions, but the quantitative pa- 
rameters were different. They suspected differences in the 
void search algorithms as reason, but we suspect that the 
effective 2-dimensional nature of the LCRS is the most likely 
cause. But their demand for using the same void search al- 
gorithm both in data and simulations seems a prerequisite 
for trustworth results. More recently, Colberg et al. (2007 in 
preparation) compared different void search algorithms and 
found that most proposed algorithm find comparable loca- 
tions and sizes of large voids. This is very likely not the case 
for the large number of small voids that fill a significant part 
of space. Therefore we shall present in this study an compa- 
rable analysis of voids both in simulations and in the data 
that explores the detailed void size distribution in depen- 
dence on the faint brightness limit of the galaxies defining 
the voids. 

W e shall use for our study the 2dFGRS (|Cole et all 
I2OO5I ) thereby coming back to the property of the self- 
similarity of the void statistics. It tells that the void size 
distribution depends on the mean galaxy separation, in such 
a way that brighter galaxies define larger voids than fainter 
ones. Even if voids in the 2dFGRS were pre viously analyzed 
l|Hovle fc VogelevI 12004 iPatiri et al.ll2006l ). this concerned 
mainly large voids and not the d etailed void statistics pro- 
posed by us previously. ICroton e t al. (2004) provided a de- 
tailed study of the void probability distribution for the 2dF- 
GRS which is related to the void size distribution but it pro- 
vides an different statistics. Essentially the void probability 
distribution is a w eighted sum over the void size distribution 
l|Otto et al.lll98d ). 

The 2dFGRS is a densely sampled survey with a com- 
pact survey geometry. This is of advantage for the question 
of the dependence of the void sizes on the galaxy magni- 
tudes defining voids. We shall derive phenomenological fits 
to the void size distribution that will be compared with sim- 
ulation results. Furthermore, we shall take up the question 
of the faint galaxies within voids. Thereby we cut both ques- 
tions, the matter content inside large underdense regions in 
the universe, and the change of galaxy properties. The color 
distributio n of galaxies in the 2 dFGRS employs SuperCOS- 
MOS data (|Hamblv et al.ll2001^ for the 7?-band. We will find 
a clear bimodality in the void galaxies s o far only studied 
in detail for the SDSS (Rojas et al. 2005; Patiri et al.ll2006l ') 
but not yet for the 2dFGRS. 

We shall evaluate our void analysis with model galaxy 
samples const r ucted from the Millenium simulation of 
ISpringel et all l|2005h and from semianalytical galaxy for- 



mation theory app lied to the numerical merger trees 
(|Croton et ahllioOa ). We analyzed specific galaxy proper- 
ties within voids and found results that can be qualitatively 
described by the model samples. A quantitative comparison 
of the galaxy color distribution and the star formation effi- 
ciency hints at certain differences between observations and 
simulations. Tentatively we connect it with specific environ- 
mental properties of galaxy formation in underdense regions. 
In particular, major mergers, galaxy harassment, tidal and 
ram-pressure stripping will not be as effective there as in 
more dense regions of t he universe (|Avila-Reese et al.ll2005l : 
iMaulbetsch et al.ll20od ). 

The outline of the paper is as follows: First we describe 
the galaxy extraction from the 2dFGRS, and in Section 3 
we provide some details of the galaxy mock data. In Section 
4 we describe our void search algorithm and in Section 5 we 
provide our results. Section 6 is devoted to a discussion and 
in the final Section we draw our conclusions. 



2 SELECTION OF 2DFGRS DATA 

We analyze t he 2d F Galaxy Redshift Survey (2dFGRS, 
ICoUess et al] (120011 )1 with about 222,000 galaxies covering 
a sky area of 1500 deg^. For the void search and statistics, 
it has the advantage of covering a sufficient large area with 
an average complet eness of 90% in redshift measurements, 
IColless et al.l (1200l|). We shall discuss the method for treat- 
ing the varying completeness below. Photometric data in 
bj st em from the APM photographic plates (Maddo x et al.l 
which represents the basis of the fiber placement for 
redshift measurements, bj and r-b and data are availab le 
through the SuperCOSMOS catalog (|Hamblv et al.l l|200lD '). 

The 2dFGRS ha s been used for i dentifying lar g e void s 
by iHovle fc VogelevI l|2004h and by IPatiri et al.1 (|2006l ). 
and for an analysi s of the void probability function by 
ICroton et all (|2004h . Our aim is to touch both points. We 
extract a large catalogue of voids, both small and large ones. 
This catalogue can be used for statistical studies of the void 
size distribution, and for identifying subpopulations within 
voids. We do not claim that all our voids are statistically 
significant. Especially the position of the large number of 
small voids are influenced by random galaxy positions. But 
a detailed look at the large voids in Figs. [T] and [2] verifies 
that both in data and simulations we catch as voids the 
same regions as the eye selects in the galaxy distribution 
(note the projection effects that influences a part of voids in 
this Figure). 

We extracted volume limited samples from the 2dFGRS 
in using the bj magnitudes since it is the color deflning the 
spectroscopic survey. Apparent magnitudes were converted 
to absolute magnitudes by using the extinction corrections 
as given in the 2dFGRS catalog for the APM bj magnitudes. 
The SuperCOSMOS mag nitudes were correct ed using the 
Schlegel extinction maps (jSchlegel et al.|[l998l ). As a proxy 
for the specifl c star formation rate we take the parameter 
as defined by 'Madgwic k et al.l l|2002l ') performing a princi- 
pal component analysis of the optical spectra. This param- 
eter also shows a strong correlation with the specific star 
formation rate obtained from the SDSS, ISol Alonso et al] 
( 20061). The paramete r 77 is correlated with the Hubble types 
( Norberg etal] |2002| ): rj < -1.4 means mainly to E/SO- 
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Figure 1. Projected distribution of the 100 largest voids in 2dF- 
GRS S3 sample for Bj < —20. Black (red in the electronic version 
of the paper) symbols are galaxies brighter than the magnitude 
threshold, and grey (orange) symbols fainter galaxies that partly 
fill the voids. Note that same bright galaxies appear inside the 
voids due to projection effects. 



galaxies, while 77 > —1.4 corresponds mostly to late type 
S galaxies. But the scatter is large and we shall interpret 
our results in terms of the morphological type with some 
care. 

All galaxies with a quality Q ^ 3 were used in order 
to have accurate redshift determinations. For determining 
th e e- and fc- c o rrecti ons we follow the method described 
in iFolkes et al.l l|l999l ) for t hose galaxies whose t ype could 
be determined, and that bv lNorberg et af] (|2002l ) for those 
which could not be determined and therefore was derived for 
a mix of galaxy types. For the Super COSMQS magnit udes 
we use the fc-correction proposed by I Cole et al] (|2005| ) for 
the r-band magnitudes. 

Eight volume limited samples where constructed both 
from coherent regions of the SGP and NGP slices of the 
2dFGRS. The SGP sample was restricted to the right as- 
cension range —2. '^10 ^ a ^ 3. '^40, and a declination range 
-33.°00 < 5 ^ -25.°50; the NGP was selected through 
10.''00 s; a «; 14.''80, -4.°00 < (5 ^ l.°00. 

The sample characteristics are listed in Table [1] The 
samples N/S 1-4 were optimized for the void search, i.e. we 
set redshift cuts that provide a large depth and a sufficient 
number of galaxies in the volume. In particular, the selec- 
tion window is chosen to be larger than typical void sizes in 
order to minimize boundary effects in the void detection al- 
gorithm. The samples have overlapping redshift ranges, but 



Figure 2. Projected distribution of 100 largest voids in a Mille- 
nium mock catalog MW3. The color coding is identical to Fig[l] 

due to different faint absolute magnitude limits -Bj,iin|3 voids 
in the different samples are defined by different galaxies. 

The samples N/S 5-6 were constructed for investigating 
the properties of galaxies inside of voids. Therefore, we have 
chosen exclusive magnitude intervals —18.4 < Bj < —19.2 
and -19.2 < Bj < -20.075 in the last four columns of Table 

m 

In order to avoid biases in the void search due to an in- 
complete sky coverage, in particular due to holes and missing 
fibers in the survey area, we take special care in getting a ho- 
mogeneously distributed completeness over the sky. In Fig. [3] 
we show the fraction of the survey area /a ( < c; ) covered by a 
completeness c below a limit c; by the dashed line. The frac- 
tion of galaxies contained in binned intervals of completeness 
c are given by the dash-dotted histogram fg{c). Both distri- 
butions follo w directly from the survey completeness masks 
provided by ICoUess etal] (|2003| ). To get a homogeneously 
sampled set, we chose a completeness cut cj, and we reduce 
all fields with c > c; randomly to the fraction cj. Then we 
get a fraction in the galaxy samples fs{ci) as 

fs{ci) = / ' /s(c)dc+ f /,(c)[l - {c-ci)]dc. (1) 

Jo J ci 

It is shown by the monotonously increasing solid line 
in Fig. O We do not care for that part of the sky covered 
only with completeness c < cj. If we chose c; = 0.6 this 

^ Note that the — 51og/i term is always included in our notation 
for absolute magnitudes Bj. 
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10609 


-18.0 


3.82 


0.59 


0.014 


0.086 


N2 


16697 


-19.0 


5.10 


2.21 


0.021 


0.135 


N3 


11749 


-20.0 


8.27 


6.64 


0.033 


0.198 


N4 


1323 


-21.0 


17.14 


6.67 


0.05 


0.2 


SI 


10343 


-18.0 


4.69 


1.07 


0.014 


0.092 


S2 


18117 


-19.0 


5.81 


3.56 


0.021 


0.139 


S3 


13964 


-20.0 


9.00 


10.20 


0.033 


0.2 


S4 


1652 


-21.0 


18.24 


10.02 


0.05 


0.2 


N5 


10389 


-18.4 


4.25 


0.80 


0.05 


0.1 


N6 


8279 


-19.2 


5.70 


1.53 


0.1 


0.14 


85 


8175 


-18.4 


5.28 


1.20 


0.05 


0.1 


86 


9516 


-19.2 


6.23 


2.30 


0.1 


0.14 



Table 1. Selection criteria for different volume limited samples 
used for void size statistics. Shown are the number of galaxies, the 
faint absolute magnitude limit Bj.iimi the mean galaxy distance 
A, the total volume V, and the minimum and maximum redshifts 
•Zmin, 2max of the volume limited samples. 




cially suited for the void search and for excluding boundary 
effects. 

The key physical processes included in the semi- 
analytical model are shock heating of gas falling into the 
dark-matter halos, gas cooling and forming of galactic disks, 
star formation in the cold phase, supernova feedback heat- 
ing of the cold phase, metal enric hment, and galaxy mergers 
leading to spheroid formation, cp. lCole et al .' (2000). As new 
feature, 'radio feedback' from AGN is implemented when a 
massive black hole forms at the centre of hot gas halo and 
suppresses further gas cooling. The radio feedback mode is 
needed in order to suppress star formation in massive dark 
matter halos. No feedback from quasar winds in included. 

The model has been shown to reproduce the observed 
2dFGRS galaxy luminosity function, the TuUy-Fisher rela- 
tion, cold gas fraction/stellar mass and cold gas metallic- 
ity/stellar mass relations for Sb/c spirals, and the global 
color magnitude relation. It also r eproduces the globa l star 
formation history of the Universe (jCroton et al.ll2()()6l ). 

The output catalog from the SAM model provides B, 
V and R band colors. We converted the B band colors to 
Bj using the relation given in Norbcrg ct al. (2002). We use 
four mock samples Ml to M4 with the same luminosity limits 
Bj,Hm = —18, —19, —20, —21 as the volume limited samples 
for the 2dFGRS. In addition we select four mock samples 
MWl to MW4 with the same magnitude limits and in ad- 
dition with the window functions of the combined samples 
N/Sl to N/S4. In addition these samples are homogeneously 
diluted to a completeness of Ci = 0.6. We do not impose the 
small incompleteness corrections for that part of the survey 
mask that is less well sampled as 0.6. If not otherwise indi- 
cated the mock samples are used in redshift space to allow 
a realistic comparison of the observed and theoretical void 
statistics. 



Figure 3. Fraction of the sky areas /a(< c) below a given com- 
pleteness c (dashed line), differential fraction of galaxies /g(c) 
sampled with completeness c (dash-dotted line) and galaxy frac- 
tion in the selected samples /g(c;) diluted to a sampling fraction 
Ci (solid line). 



concerns only about 5% of the area. On the other hand, we 
have about 75% of galaxies in the samples fs{ci). This is a 
straight forward procedure that is adopted to our specific 
problem. 



3 SIMULATED GALAXIES 

For a quantitative comparison with CDM models of galaxy 
formation, we compare the voids statistics and galaxy void 
properties with tho se predicted by the s emianalytical galaxy 
formation model of ICroton et al. (|2006l ) . 

This model uses the dark matter halo distribution and 
mass accretion histor ies found in the Millenium simulation 
l|Springel et al.l [20051 ). and applies a semi-analytic formal- 
ism to describe the evolution of galaxies in these halos. 
The Millenium simulation covers a box of sidelength 500 

Mpc with 2160^ particles, run in a concordance cosmol- 
ogy {ilm = 0.25 and as — 0.9). The large volume is espe- 



4 FINDING VOIDS 

Although it is generally agreed upon that voids constitute 
underdense regions in the galaxy distribution, the precise 
definitions of voids varies considerably in the literature. 
Sometimes a threshold in the cosmic density field is im- 
posed to define voids, either from the dark matter density, 
or constructing a continuous density field from the galaxy 
distribution. In either case, we must choose a smoothing 
scale to avoid strongly varying densities and corresponding 
irregular voids. A variant is the de finition of wall galax - 
ies from the local density field, cp. lEl-Ad fc PiranI (|l997D 
and lHovle fc VogelevI (|2002l '). Others define voi ds by regions 
empt y of galaxies over a certain magnitude (|Muller et al.l 
I2OOOI: [ Patiri et al. IJ2OO6) or of halos over a mass limit 
l|Gottl6ber et al.l l2003l. 

In many previous void searches, the geometry of the 
voids is specified to be spherical, or a union of overlapping 
spheres. Although from gravitational instability theory it 
follows that voi ds become more spherical during evolution 
as they expand (Ilckelll984l : lvan de Weveaert fc van KampenI 
1 1993) , in real data and during intermediate evolution stages, 
they can have a wide range of geometries (|Shandarin et al] 
l2004h . 

Voids can be seen as underdense perturbations in the 
initial Gaussian dark matter density field that expand and 
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become non- linear. Using the top-hat spher i cal co llapse 
(expansion) model, ISheth fc van de Wevgaer'^ (|2004l ) have 
shown that at the shell cr ossing time, the linear mean den - 
sity should he 5 = —2.41 jSheth fc van ^ Wevgaertll2004h . 
In a recent paper lFurlanetto fc Piran ( 20061 ) have shown how 
this can be related to the mean galaxy density, using an 
Halo Occupation Distribution (HOD) model. Although this 
definition is appealing, the corresponding galaxy density de- 
pends on the size of the smoothing length. It is therefore not 
obvious so far how to relate uniquely the galaxy distribution 
to the underlying dark matter distribution. 

Therefore in this paper we define the voids as regions 
being empty of objects larger than some given absolute mag- 
nitude Bj. As the more massive and brighter galaxies will 
be found in filaments and clusters, we use them as tracers of 
the cosmic web. To investigate the self-similarity of the void 
distribution, we analyze a sequence of increasingly bright 
galaxies B.j = -18,-19,-20,-21, the samples N/Sl-4 in 
Table m 

W e used the void finder described in (|Mtiller et al.l 
|2000[) . It is a grid-based void finder that locates the largest 
empty cube on the density grid. Next it looks for neighbor- 
ing empty density layers along all six sides of the cube. The 
largest compact empty layer consisting of a base square layer 
and extensions along the four square boundaries is added to 
the void if its volume is a factor / = 0.67 larger than the 
previous boundary layer of the cube. This process is iter- 
ated along all six boundaries of the base cube, each time 
requiring that the additional layer is a boundary of the pre- 
vious void boundary and exceeds this by the given factor 
/. This factor is also imposed in extending the base square 
layer along their for sides. This procedure provides almost 
convex empty voids. The value / is the only free parame- 
ter in our void search algorithm. Its value is chosen so as 
to avoid narrow tunnels that go out of the base voids into 
the galaxy distribution. For the present analysis, we provide 
only results for the base square voids. We have analysed both 
2dF-data and simulations also including the void search with 
extensions, with gives similar void size distributions, scaling 
relations and an reproduction of the data by the simulations. 
The reason for restricting on the base voids lies in its use for 
looking for faint galaxies laying in central regions of large 
voids. 

For the analysis, a grid of resolution of 1 Mpc was 
chosen. The void finder is useful in that it is faster than 
looking for the largest empty spheres in the galaxy distribu- 
tion. Although voids do not tend to be cubical, the extension 
along the boundaries and a sufficiently small grid allow for 
realistic voi d geometr ies. Comparison with the spherical void 
finder of iPatiri et al. I (2006 ) shows, that for large voids we 
are mainly interested in, the void finder selects the identical 
void centers and sizes. 

In our analysis we derive the size distribution of voids 
up to small values, thereby covering in most cases over 90 % 
of the space with voids. For the comparison with other void 
finder, we use and effective radius R = {SV/iny^^ . As statis- 
tics we use the volume weighted void size distribution 
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f{R)dR, 



(2) 



vantage of allowing a robust void size distribution over 2 
to 3 order of magnitude in F{> R). The strongly varying 
volumes of large voids then appear at the low abundance 
end of the distribution. The more common void probability 
function Po(^) is given by a weighted sum o ver the dif- 
ferent ia l size distributi o n fjR ) — —dF/dR, cp. lOtto et al.l 
d 19861 ): iBetancort-Riid (|l990l ). The employed void size dis- 
tribution F{> R) does not suppose a pre-specified geometry 
of the voids, and it determines the maximal empty region in 
the galaxy distribution, starting from large voids and going 
to smaller ones. Thus it seems to be most sensitive to the 
large-scale distribution of galaxies in the cosmic web. 



5 RESULTS 

5.1 The void size distribution 

First we investigate the void sizes found in the observed 
volume limited samples. We found maximum voids with ef- 
fective radii of 24 Mpc and 22 Mpc base sizes in the 
sample N4 and S4, resp. There are 156, 70, and 2 voids over 
12 Mpc base length in samples N4, 3, 2, and 1, resp.; 
and 199, 116, and 2 voids in S4, 3, and 2, resp. (no such 
large void are found in Nl or SI). Voids over 6 h-^ Mpc ef- 
fective radius are much more abundant, there are 916, 1250, 
321, and 66 voids in samples N4, 3, 2, and 1, resp.; and 1364, 
1786, 538, and 129 voids in samples S4, 3, 2, and 1, resp. Fig. 
[T]shows the 100 largest voids as found in the sample S3. For 
comparison Fig. [2]shows the same for a Millenium mock cat- 
alog with the same magnitude range and an identical survey 
mask (i.e. MW3, see below). Obviously, both 2dFGRS- and 
mock samples look very similar, and they show a compara- 
ble distribution of large voids. Large voids are more abun- 
dant at larger distances from the observer. This is an effect 
of the survey window that strongly restricts large nearby 
voids. Besides this geometric effect, large voids are evenly 
distributed in data and model distributions, and void sizes 
are comparable in both samples. 

More quantitative results can be seen in the cumulative 
void size distribution shown in Fig. |4] There we combine 
the S/N samples from the 2dFGRS to get better statistics. 
If taken separately we get similar curves for the north and 
south galactic pole regions. We take the difference for the 
separate distributions as a measure of the uncertainty in 
the samples shown by the error bars on the binned sample. 
Most obvious is the size dependence of voids on the galaxy 
samples, the largest voids among bright galaxies with Bj > 
20, 21 have an effective radius over 15 Mpc//i. The statistics 
of such large voids is still quite restricted, but smaller voids 
lead to size distribution with small error bars. The general 
form of the distribution is self-similar. We fit it with modified 
exponential distributions 



F{> R) = exp 



\si\) ^{S2\ 



(3) 



defined as the cumulative volume of the survey covered by 
voids of radius larger than R. This definition has the ad- 



with 4 parameters, two length factors si and S2, and two 
powers pi and P2- The fits are by ^^-minimization and al- 
ways provide excellent representations of the observational 
data covering more than 2 orders of magnitude in -F(> R). 
The main scaling of the galaxy samples is described by the 
mean separation between galaxies, the scale A. This value 
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and all fit parameters are given in Table (2] The first ex- 
ponential describes the statistics of small and intermediate 
sized voids. This part is typically described with an linear 
fall-off in the exponential function, i.e. pi ~ 1. For describ- 
ing the cutoff of the void size distribution for large void 
radii, we need the second factor in the exponential distribu- 
tion with a higher power of the void radius, typically with 
P2 ~ 3, i.e. the abundance of large voids is cut off with the 
void volume entering the exponential function. The length 
factors si and S2 are both of order unity. After multiplica- 
tion with the mean galaxy separation they provide the void 
radii where the two exponential laws dominate the size dis- 
tribution. There are no ways to fit the size distribution with 
a single exponential distribution. 

In Fig. [S]we show the normalized void size distribution 
of the four 2dFGRS data sets, i.e. we divide the void sizes R 
by the mean separation between galaxies A. Obviously, the 
void size distributions are self-similar in this variable. The 
scaled void size distributions coincide at radii R < 1.5A. 
This scaling is similar to the dependence o f the correlation 
length of g alaxies and groups as found by iBahcall fc WestI 
l|l992h and I Yang |2005f ). Note that the samples N/S4 

are strongly restricted by the survey geometry. The mean in- 
tergalaxy separation of about 18 Mpc is comparable to 
the thickness of the 2dFGRS slices of about 45 Mpc at 
the far side of the selected volume, even if the length of the 
slices much exceeds the void sizes. Effectively the statistics 
of these voids, shown by the dash-dotted lines in Fig. (5] is 
cut off at an volume filling fraction of 2 %. Most interesting 
in the scaled void size distribution is that the exponential 
cutoff in the void size distribution depends on the galaxy 
sample defining the voids. Voids among fainter galaxies have 
a larger tail in the distribution of void radius over mean 
galaxy separation, R/X. In this variable, the largest voids 
are among the faint galaxy samples N/Pl and 2. Gener- 
ally the void radii R/X extend over larger ranges than voids 
in Poisson samples. The voids are larger than in a random 
galaxy distribution, and the size distribution shows a sharp 
cut-off what is an indication the the voids are hitting the 
galaxy walls in the large cosmic web. 

We compare the void size distribution of the 2dFGRS 
with the mock samples from the Millenium simulation in 
Fig. m Here we show the samples N/S2 with the diamond 
symbols and error bars. The solid histograms and the fit 
from Table [2] shows the results for the Millenium galaxy 
mock sample restricted to the 2dFGRS survey window (the 
sample MW2) with Bj < —19. In addition we show the 
void statistics for the complete Millenium volume (sample 
M2) in real (dotted line) and redshift space (dashed line). 
These two distributions are quite similar, but systematically 
shifted to larger void radii by about 5% larger in redshift 
than in real space. This is an effect of the overall expan- 
sion of the underdense regions in comoving coordinates, i.e. 
forground and background void boundaries show systematic 
deviations from the Hubble flow. The mock sample MW2 
in the 2dFGRS-window is shown in redshift space. Within 
the error bares, it well reproduces the void statistics of the 
N/S2 sample. There is a slight tendency of overpredicting 
the abundance of large voids in the simulations. It should 
be noted that the void sizes R in the 2dFGRS survey win- 
dow are nearly 20% smaller than in the 500 Mpc . Still 
for the 2dFGRS window geometry, finite volume effects are 
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Figure 4. Cumulative volume weighted void size distribution 
F(> R) in the 2dFGRS samples N/Sl (triangles), 2 (diamonds), 
3 (squares) and 4 (stars). The analytic fit according to eq. (3) 
provide excellent descriptions of the data. 
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Figure 5. Normalized void size distribution F{> D/X) as func- 
tion of the void diameter divided by the galaxy separation A in 
the 2dFGRS with same line symbols for samples N/Sl-4 as in 
Fig.H 

important for all the range of the void statistics. The general 
agreement of the model and 2dFGRS void statistics holds 
true for the other data N/Sl, 3, and 4; and the mock sam- 
ples MWl, 3, 4, as the parameters of the flt curves in Table 
[2] demonstrate. 

Most interesting is the dependence of the void size dis- 
tribution on the limiting magnitude, shown in Fig. [7] for the 
full statistics of the Millenium box. This is the complement 
to Fig. [5] Below void radii R < 1.5A, the void size distribu- 
tions coincide for different samples. This means void sizes 
depend strongly on the mean galaxy separation, and for a 
long range they are proportional to A. Beyond this radius, 
R > 1.5A, void sizes depend strongly on the magnitude cut. 
Voids among faint galaxies have an more than three times 
larger radius R than the mean intergalaxy separation, see 
the results for Bj > — 18 galaxies (dotted histogram). In the 
contrary, voids among bright galaxies with Bj > —21 have 
radii only up to 2A. This is an expression of lPeeblesI (|200ll ) 
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catalog 


A 


Pi 


P2 




S2 




Mpc/h 










N/Sl 


4.26 


0.8 


3.0 


1.02 


1.49 


N/S2 


5.45 


1.1 


4.5 


0.86 


1.80 


N/S3 


8.64 


1.0 


4.5 


0.93 


1.47 


N/S4 


17.69 


0.6 


4.0 


0.67 


0.86 


MWl 


4.16 


0.7 


2.0 


1.70 


1.26 


MW2 


5.55 


0.8 


3.0 


1.10 


1.50 


MW3 


9.13 


1.0 


4.0 


0.86 


1.46 


MW4 


22.21 


0.5 


4.0 


0.43 


0.70 


Ml 


3.58 


0.9 


3.0 


1.36 


2.07 


M2 


4.83 


1.0 


3.5 


1.23 


2.02 


M3 


8.02 


1.1 


4.0 


1.18 


1.71 


M4 


19.60 


0.6 


3.5 


3.79 


1.17 



Table 2. Fit parameters of void statistics for volume limited 
samples of the 2dFGRS and of galaxy samples (NGP and SGP) 
and in the Milleninum simulation (MW with window selection 
from 2dFGRS, M for the complete box). The different samples 
are characterized by the mean galaxy distance A. 




Figure 6. Comparison of the S/N2 void size distribution of the 
2dFGRS with the data from the Millenium mock catalogue MW2 
modelling the selection window of the observations (thick his- 
togram and solid fit curve). The dashed and dotted lines are the 
void size distributions of the complete Millenium box in redshift 
and real space, resp. 



void phenomenon: There are relatively large voids among 
faint galaxies. And otherwise, voids among bright galaxies 
are also underdense in faint galaxies which are found pre- 
dominantly at in the outer regions of the voids from bright 
galaxies. 

The self-similarity of the void size distribution that is 
seen in the proportionality of the void size distribution to 
the mean intergalaxy separation becomes also obvious in the 
radii of voids that cover 25 %, 50 %, and 75 % of the space, 
i.e. the median and the quartiles of the void size distribution 
of Fig. m The filled squares in Fig. |5]show these percentiles 
of the Nl-4 and Sl-4 volume limited samples of the 2dF- 
GRS. The differences in the mean galaxy separations and 
void percentiles between neighboring data points show the 
uncertainty due to cosmic variance. They are of the same or- 
der as the symbol sizes. These percentiles are well described 
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Figure 7. Normalized void size distribution in the complete Mil- 
lenium box for magnitude limits Bj < —18, —19, —20, —21, i.e. 
the samples Ml-4 from right to the left. 



by linear fits to the lower 25% quartile i?25, the median i?5 
, and the upper 75% quartile f?75 , 



Ro + vx A, 



(4) 



where an initial value Ro = 1.1±.2, 1.8±.4, 3.5±.2 Mpc 
and aslope v = 0.23± .02, 0.51 ± .04, 0.79±.01 for the 25%-, 
medium and 75%-percentiles are found by linear regression 
of the data. The squares correspond to the mock samples 
MWl-4 from the Millenium galaxy catalogues with limiting 
magnitudes of B < —18, —19, —20, and— 21. For the 3 fainter 
samples, they fall completely on the observed similarity rela- 
tions. The triangles are the samples from the full Millenium 
box Ml-4. For samples with A < lQh~^ Mpc they lie about 
10 percent above the data and simulation data in the 2dF- 
GRS window. These samples are not infiuenced by boundary 
effects. The large voids for A > 15/i~^ Mpc show even larger 
differences between the samples in the simulation box and 
in the window volume. This similarity r elations were previ - 
ously found in the LCRS void ana lysis ('M iiller et alll2000l ) 
and in CDM simulations l|Benson et al. .200j ) but there with 
less statistics, i.e. larger uncertainties. 



5.2 Properties of void galaxies 

Now we describe our results concerning the void galaxies. 
We ask whether galaxies with magnitudes fainter than the 
lower absolute magnitude limit -Bj^nm of the volume lim- 
ited samples for the void search are populating the voids 
and whether their properties differ from average galaxies. 
Because by definition the voids are sparsely populated we 
stack all voids over a minimum void size Rv ^ Rvum and 
look for all void galaxies within a scaled radius /_b = R/Rv- 
We investigate voids identified in the volume limited sam- 
ples N/S5 and 6 with magnitude limits B,j G [-19.2, -18.4] 
and Bj € [-20.075, -19.2]. We take all voids larger than 
the limiting diameter Rvum = Wh~^ Mpc for the analysis. 
To attain better number statistics we include all galaxies in 
the detected voids, also those that had been removed while 
correcting for nonuniform completeness during the analysis 
of the void size statistics. In this way we obtain a sample of 
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A (li' Mpc) 

Figure 8. Median and quartiles of the volume covered by voids 
as function of the mean intergalaxy separation. The full dots stem 
from 2dfGRS data, the squares from mock samples MWl-4, and 
the triangles from the complete Millenium box samples Ml-4. The 
lower 25% percentile, the medium and the upper 75% percentile 
are described by linear fits, curves from below. 



661 galaxies in the fainter and 1222 galaxies in the bright 
sample. Although the luminosity function is a decreasing 
function with increasing luminosity, we obtain more galax- 
ies in the brighter sample, because their brightness allows a 
deeper redshift range and a larger volume (cp. Table [1]). We 
choose to analyze separately the central parts of voids with 
fn = 0.3 and the overall properties of void galaxies with 
fR = 0.6. 

In Fig. [5] we show the normalized color distributions of 
void galaxies within /r = 0.6 as histograms with Poisson 
errors together with double Gaussian fits obtained with y^- 
minimization as solid histograms and lines, respectively. As 
control sample for average galaxies, we randomly choose an 
equal number of spherical regions in the 2dFGRS volume as 
the number of voids. The number of galaxies is given in the 
upper part of the diagrams. We found a clear bimodality 
both in the void galaxies and in the control sample for both 
brightness bins (taken NOP and SGP samples together). 
This is the first time that this color bimodality has been 
shown for void galaxies in the 2dFGRS. Both for the brighter 
(lower panel) and fainter (upper panel) samples, inside voids 
the fraction of blue galaxies increases on account of the frac- 
tion of red galaxies. In as far as the color is a measure of the 
morphology of the galaxies, the changed proportion of blue 
and red galaxies is a continuation o f the morpholog y density 
relationshi p of cluster galax ies by iDressleil (|l980t ) and dis- 
cussed by iHogg et al.l l|2004 ) for SDSS galaxies estimating 
the density on a local scale of 1 Mpc projected on the 
sky and by 1000 km/s in radial direction. Fig. |9] shows that 
the suppression of the red sequence inside voids is stronger 
for the brighter galaxies. The ratio of void galaxies in the 
red and blue clouds are almost independently of the bright- 
ness, while the field population has a stronger red sequence 
population in the brighter bin, i.e. there is a stronger effect 
of star formation strangulation than in voids. 

For the brighter galaxies, the color distribution of the 
void sample has been analyzed for difi'erent fractions of the 
void size, /_r = 0.3,0.6 (Fig. I10|) . It is interesting that the 
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Figure 9. The magnitude dependence of the (Bj — R) color 
distribution of galaxies in void centers R/Rv ^ 0.6 (solid his- 
togram with 1(7 Poisson error bars) and double Gaussian fit (solid 
line). Top: The fainter and bottom the brighter magnitude ranges 
given in the upper right corner of both panels). The dash-dotted 
lines are the control samples in the complete survey volume. The 
galaxy numbers are shown in the upper left corners. 



blue cloud seems to show a blue shift by about 0.1 mag, 
while the red sequence remains fixed at about B — R = 1.2. 
So there is a weak tendency that the blue galaxies in most 
underdense regions are not only more abundant, but also 
tend to be in average bluer and maybe younger than galaxies 
in the same magnitude range in the control sample. The 
faint bin does not contain enough galaxies to plot the color 
distribution for more centrally confined galaxies. 

For the same selection criteria we show the distribution 
of the rj parameter as a proxy for the specific star formation 
rate in Fig llll For both the bright and less bright sample, 
the star formation rate of void galaxies is enhanced, although 
this is only marginally significant for the lower magnitude 
bin. 

From the semi-analytical model in the Millenium simu- 
lation we get qualitatively similar color distributions of the 
void galaxies. We show in Fig. [12] the [Bj — R) color dis- 
tribution of the faint (above) and bright galaxies (below). 
Again the distributions can be well fitted with double Gaus- 
sian distributions. In the model samples, the mean colors of 
both blue and red galaxies are the same for voids and the 
field, i.e. the void environment has an infiuence on the rel- 
ative abundance of the both galaxy species, but not on its 
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Figure 10. The radial dependence of tlie (Bj — R) color distribu- 
tion of the brighter void samples. The solid distributions show the 
R/ Rv ^ 0.6 and the dashed distributions the central R/Rv ^ 0.3 
void regions. The dash-dotted line repeats the comparison sam- 
ple, the galaxy numbers are again given in the upper left corner. 
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Figure 12. The magnitude dependence of the color distribution 
{B j — R) of void galaxies in the Millenium mock galaxy sample. 
The solid distributions show the void and the dash-dotted distri- 
butions the control samples. Top: The faint and bottom the faint 
galaxy sample. 
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Figure 11. Top: Distribution of the spectral parameter r] in 
voids (R/Rv ^ 0.6) (solid lines) and in void centers R/Rv ^ 0.3 
(dashed lines), the control sample is given by the dash-dotted 
line. Top: The faint and bottom the bright galaxy sample. 



colors. Comparing Fig. [9] and Fig. I12l we denote some quan- 
titative diflterences. The suppression of the red peak inside 
voids is more pronounced for the faint magnitude bin, and 
altogether the red galaxies have a much tighter color range. 
The peak of the red sequence of the brighter galaxies is much 
sharper than in the observed sample, and it lies at a slightly 
less red color, Bj — R « 1.1. So although the location and the 
width of the blue cloud is in good agreement with the data, 
the semianalytical models have too many red galaxies with 
a too sharp color distribution. Furthermore, the magnitude 
dependence of the ratio of red and blue galaxies shows sig- 
nificant differences to the data. The faint galaxies show the 
stronger suppression of the red sequence than the brighter 
magnitude bin. A similar discrepancy of the fraction of red 
galaxies inside galaxy groups of the SDSS redshift survey 
against the Millenium si mulation semianalytic ga laxy cata- 
logues has been found bv lWeinmann et all (|2006l ). 

Finally, we compare the radial distribution of the frac- 
tion of late (S) and and early type (E) void galaxies in the 
2dFGRS and in the Millenium simulation in Fig. [I3l In the 
2dFGRS we employ the spectral parameter rj for roughly 
distinguishing early and late galaxy types, with rj ^ —1.4 
for E and rj > —1.4 for S galaxies. In the semianalyti- 
cal model, S galaxies are defined by a mean star forma- 
tion rate logSFR/Msun/yr > 0, in the contrary E galaxies 
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Figure 13. Radial distribution of the fraction /typo of early (E, 
lower spherical symbols) and and late (S, upper triangle symbols) 
in dependence of the scaled void radius. Solid symbols stem from 
the 2dFGRS and open symbols from the Millenium mock galaxy 
sample both in the redshift range —18.4 > Bj > —20.075. 



by log SFii/A/sun/yr ^ 0. These definitions lead to a frac- 
tion / = 0.52 of E-type galaxies botii in the 2dFGRS and 
the Millenium catalogue. We show the fraction of E- and 
S-galaxies within voids of radius larger then 13 Mpc 
and for the combined samples N/S 5 and 6, i.e. the range 
-18.4 < Bj < -20.075, solid symbols for the 2dFGRS and 
open symbols for the Millenium galaxies. Obviously both 
distributions are quite similar. They coincide within the la 
error bars at mean radii values but the 2dFGRS shows a 
stronger radial variation of the spiral and elliptical fraction 
There is the general tendency that the fraction of E-galaxies 
in observed voids is smaller than in the simulations, this 
tendency extends over the total radial range of voids. Fur- 
thermore, in the 2dFGRS, the difference in the E/S-fraction 
extends beyond the void volume, i.e. the differences in the 
star formation activity is more pronounced and extends over 
larger scales in the data than in the comparison semianalyt- 
ical model. 



6 DISCUSSION 



ing into account the transition to the 3-dimensional survey 
geometry. 

The scaling relation (i.e. the constant of prop ortional- 
ity) of void sizes i? oc A is flatter than found bv , Benson et al.l 
(|2003t ). We suspect this is due to our procedure of identifying 
both large and small voids with the same algorithm, not just 
looking inside of wall galaxies of fixed local overdensity as 
other void finders do. The scaling relation is a representation 
of the self-similarity of the void phenomenon. In distinction 
to earlier analyses, we have a large sample of voids and there- 
fore a good statistics in determining the void size distribu- 
tion. The flat scaling relation is a consequence of relatively 
small voids among the bright galaxy samples with large A 
and conversely, comparatively large voids among faint galax- 
ies. This is a new property well established in the densely 
sampled 2dFGRS and in the large Millenium simulation. It 
is the sa me void p henom enon discussed in the stimulating 
paper by iPeeblesI ( 200 ll ). However, he claims a significant 
discrepancy of the void phenomenon and the standard hier- 
archical clustering phenomenon. Contrary to this claim we 
established void statistics that are well reproduced by the 
Millenium galaxy catalogues. As shown in Fig.[7]the size dis- 
tribution of voids in simulated galaxies is well described by 
a self-similar distribution for small and intermediate sized 
voids, and a more extended tail of large relative void sizes 
with respect to the mean galaxy sep aration. This m eans that 
the void phenomenon discussed bv IPeeblesI (|200lh is a gen- 
uine consequence of the underlying LCDM model. An earlier 
demonstration of the void phenomenon within the hierarchi- 
cal CDM cl ustering model by quite different methods was 
provided bv'Mathi s fc White! (|2002l ). 

The new results concern the precise verification of the 
void scaling relations for the 2dFGRS survey and the de- 
tailed fitting of the void size distribution by a two-branch 
exponential distribution. Both properties of the 2dFGRS are 
surprisingly well reproduced in the semianal ytical galaxy 
form ation scheme of the Millenium simulation ( Croton et al.l 
The comparison is done in the 2dFGRS survey mask 
and in redshift space. This is essential for the precision com- 
parison since we found coinciding void size distributions 
both for small and median void sizes, but an increase of the 
largest voids in redshift space by about 10%, and an overall 
reduction of the void size distribution in the 2dFGRS survey 
mask as compared to the complete 500 Mpc simulation 
box. 



6.1 Void size distribution 

We found a clear demonstration of the similarity relation of 
void sizes in the 2dFGRS, iZvoid oc A, first est ablished em- 
pirica lly both in observations and in models l|Miiller et al.l 
l200d ). There the much smaller LCRS was analyzed and the 
analysis was done in 2-dimensional slices due to the geome- 
try of this survey. Later similar relations were found for voids 
in the da rk matter density fi eld using the void probability 
function JSchmidt et al ]|200ll). among CDM-halo catalogues 
l|Arbabi- BidgoU fc Mxillcr"2002'). and among semianalytical 
galaxy catalogues (Benson ct al. 2003). Here we could pro- 
vide an analytical description of the void size distribution 
given by Eq. (3). The similarity relation Eq. (4) is a direct 
consequence of such a void size distribution. The scaling of 
the median sizes is comparable to our earlier findings tak- 



6.2 Void galaxies 

We established a clear representation of the color bimodality 
of void galaxies and an larger fraction of blue galaxies in the 
2dFGRS voids. An increase of the fractio n of blue galaxies 
in voids has been previously esta blished bv l Hovle fc VogelevI 
(|2004l ) for the 2 dFGRS a nd by R qias et al.l (|2004 12005 ^) for 
the SDSS. Also. lCroton'et al., (,2001 )" found a dependence of 
the blue fraction on t he large-scale den s ity en vironment of 
a scale of 8/i"^ Mpc. iKauffmann et all {2004 ) investigated 
the SDSS for the environment dependence of star formation 
activity of galaxies on a more local scale. They found that 
the most sensitive environmental dependence of galaxies is 
the SFR, which is strongest for galaxies with stellar masses 
< 3 X 10^° Mq, but no much environmental dependence on 



Voids and Void Galaxies 11 



the large scale density was found when the small scale den- 
sity was specified. 

A bimodality of the color distribu tion of void g alaxie s 
in the SDSS has been established by iPatiri et alj (|2006l ). 
Our results for the 2dFGRS are new and extend these ear- 
lier findings. As in that paper, we confirm the stability of 
the red sequence in void regions, but we find a slight blue 
shift of the blue cloud. Furthermore, there is a weak ten- 
dency of void galaxies to have a higher SFR, and for more 
early type (E) than late type (S) galaxies. A reason for the 
bl ue slight shif t of th e blue cloud established here but not 
bv lPatiri et al.l (|2006l ) in the SDSS analysis may be that we 
analyzed the central void regions, while they considered the 
whole void volume. Also the 2dFGRS is denser sampled, and 
therefore we have a better statistics of void galaxies. While 
the SDSS colors have a b etter photometric accuracy, it was 
shown by iNorberg et al. I 11)02) that there is only a small 
non-linearity between the APM and SDSS colors. Some dif- 
ference might be produced by differences in the adopted k- 
and e-corrections. 

Color distributions at higher redshifts show a similar bi- 
modal behavior as our results. By comparing the SDSS color 
distribution with DEEP2 results. iBlantonI (|2006| ) found that 
the blue cloud is shifted bluewards by 0.6 mag at redshift 
z — 1. Of course, this shift is stronger than what we find 
here, but it indicates a similar trend. Therefore a natural 
interpretation would be a slower evolution of structures in- 
side voids and a younger galaxy population as compared to 
the field. 

If the color shift might be due to a younger population 
in underdense regions, we would expect the shift in the blue 
peak to be stronger for the lower magnitude bin. However 
due to the low number of void galaxies in the fainter mag- 
nitude bin, we cannot establish such a tendency. 

The color shift is of similar amount than the color un- 
certainty of 5m — 0.09. But an uncertainty would basically 
increase the broadness of the blue cloud galaxy distribution, 
and not leading to a systematic shift. Therefore we con- 
clude that the precise color distribution of blue void galaxies 
should be further investigated. 



7 CONCLUSIONS 

In this paper we have studied properties of voids in the 
2dFGRS and compared them with those predicted by semi- 
analytical models of galaxy formation. In particular, we were 
interested in the distribution of void sizes in the galaxy dis- 
tribution as a function of the faint limiting magnitude of 
the sample. We established an almost linear dependence of 
the void size distribution on the mean separation A of the 
galaxies in the sample. A similar relation was found for the 
magnitude dependence of the correlation leng t h of galax- 
ies a nd galaxy groups (|Bahcall fc WestI 1 19921 : lYang et al] 
I2OO5I ). It is a natural consequence of the halo model of 
gravitational clustering and the resulting void statistics 
ijTinker. Weinberg fc Warrenll2006l ). 

In addition to the self-similarity, we have found that 
the voids among faint galaxies extend to relatively larger 
scales when devided by the mean void size. In the scaled 
radius R/X, the voids size distribution of faint galaxies has 
a longer tail at large radii than those for brighter galax- 



ies. This indicates that the faint galaxies trace the general 
spatial distribution of the cosmic web of brighter galaxies. 

For the 2dFGRS, we find more large voids at larger dis- 
tances from the observer. This is due to the conelike struc- 
ture of the survey. In our mock sample, we use exactly the 
same geometry of the observed samples. Thereby we found 
that the number density of the largest voids can be under- 
estimated by up to 20%. 



We confirm the prev i ous findings by Grogin fc Geller 
(2000); H ogg et all (l2004lj : lRoias et all l|2005l ): ICroton et al 



(2005i '): iPatiri et al.1 (|2006l ). that in lower-density regions, the 
galaxy population is dominated by blue actively star form- 
ing galaxies. We fitted the Bj — R color distribution of 2dF- 
GRS void galaxies by double Gaussian distributions. There 
are significantly more void galaxies in the blue cloud than 
in the general field, and the red sequence is strongly sup- 
pressed. In addition, we have found the indication that the 
blue population of galaxies in void centers is slightly bluer 
than the field population, not just more numerous. The shift 
is only minor and alm ost of the same or der as the maximal 
quoted uncertainty bv lCole et al.l (|2005l ). so it remains un- 
clear how robust is this effect. It is not reproduced by the 
semi-analytical models for galaxy formation of the Mille- 
nium simulation. Also the red sequence in the Millenium 
catalogue has a significantly tighter spread. This may due 
to a too sharp cutoff of the star formation activity in the 
semianalytical models which seems to be unrealistic. 

The radial distribution of the ratio of early and late 
type galaxies inside voids coincides quantitatively for the 
2dFGRS and the Millenium catalogue, but in the data, the 
fraction of E-galaxies stays smaller (and the fraction of S- 
galaxies larger) further outwards beyond the boundary of 
the voids. The abundance of star-forming galaxies is higher 
in the 2dFGRS voids than in the semianalytical model. 
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